SF  298  MASTER  COPY 


KEEP  THIS  COPY  FOR  REPRODUCTION  PURPOSES 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  NO.  0704-0188 


^am«rinn  ann  ma  informiiion  <s  esiinitted  10  avarag*  i  nour  o«r  rasoonsa,  tnciudmg  ina  tima  for  reviewing  mstructionj,  saarcning  axisimg  data  iourca*, 

‘-'iiaction  nt  tM»  data  ne«oad.  and  completing  and  reviewing  me  collection  of  information  Send  comment  regarding  this  Durden  estimates  or  any  other  aspect  of  this 

^avis  Hionwav%  Washington  Heaoduaners  Services,  Oirectoraie  for  information  Ooerations  and  Repons.  1215  Jefferson 

^  *  Q  ay  iJuiie  i204  Arlington.  VA  22202-4302.  ano  to  tne  Office  of  Management  ana  Budget.  Paoerworx  ReducNcn  Proiect  10704.0188),  Washmoton.  DC  20503 


1.  AGENCY  USE  ONLY  (Leave  PlanK} 


2.  REPORT  DATE 

December  1996 


ana  Budget.  Paoerworx  ReducNcn  Proiect  10704.0188),  Washington.  DC  20503, 

3.  REPORT  TYPE  ANO  DATES  COVERED 
'  Technical 


4,  TITLE  ANO  SUBTITLE 


5.  FUNDING  NUMBERS 


nergy  Absorption  for  Short  Duration  Impacts 


DAAH04-96- 1-0082 


6.  AUTHOR(S)  - 

M.  Ziejewski,  B.  Anderson,  M.B.  Rao  and  M.  Hussain 


7  PERFORMING  ORGANIZATION  NAMES(S)  ANO  ADORESS(ES) 
Center  for  Multivariate  Analysis 
il7  Thomas  Building 
Department  of  Statistics 
Penn  State  University 
Jniversity  Park,  PA  16802 


9  SPONSORING  <  MONITORING  AGENCY  NAMEiSl  AND  ADDRESSlESl 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 


96-14 


10.  SPONSORING /MONITORING 
AGENCY  REPORT  NUMBER 


I  .S,  ,Armv  Research  Office 
P.O.  Bo,x  l2211 

Research  Triangle  Park,  NC  27709-221 1 


11.  SUPPLEMENTARY  NOTES 

The  Views,  opinions  and/or  findings  contained  in  this  report  are 
an  official  Department  of  the  Army  position,  policy  ortfecision. 


those  of  the  author! s)  and  should  not  be  construed  as 
unless  so  designated  by  other  documentation. 


12a.  DISTRIBUTION  /  AVAILABILITY  STATEMENT 


12  b  Distribution  CODE 


.Approved  for  public  relea.se;  distribution  unlimited. 


19970212  093 


13.  ABSTRACT  (Max/num  ^00  wor(Js; 


The  objective  of  this  study  is  to  relate  energy  absorption  characteristics  to  selected  ma¬ 
terial  properties  and  to  establish  a  methodology  that  allows  one  to  determine  some  of  the 
material  properties  for  meiximum  energy  absorption.  The  finite  element  program  DYNA-3D 
and  its  ztssociated  pre  and  post  processors  were  used.  The  model  used  is  a  hollow  square 
colunm.  Five  properties  of  the  materials  were  included  in  the  analysis:  (i)  Density  (ii) 
Elastic  Modulus  (ill)  Tangent  Modulus  (iv)  Yield  Strength,  and  (v)  Poisson  Ratio.  The  Re¬ 
sponse  Surface  Method  in  conjunction  with  the  canonical  analysis  were  employed  to  locate 
the  optimum  or  near  optimum  levels  of  the  properties  and  then  to  determine  the  equation 
of  the  response  surface  in  an  area  near  the  vector  of  optimum  levels.  For  the  given  levels  of 
three  out  of  five  material  properties  used  in  the  study,  one  can  calculate  the  remaining  two 
material  property  levels  to  achieve  the  near-optimad  energy  absorption. 


15.  NUMBER  IF  PAGES 
7 


16.  PRICE  CODE 


SECURITY  CLASSIFICATION  19.  SECURITY  CUSSIFICATION 
OR  REPORT  OF  THIS  PAGE  OF  ABSTRACT 

UNCLASSIFIED  UNCLASSIFIED  lJNri.A.<:.«;rFTFn 


20.  LIMITATION  OF  ABSTRACT 


NSN  7540-01 -280*5500 


xr  msPECTED  3 


UNCLASSIFIED 


Standard  Form  298  (Rav.  2-89) 
Prescribed  by  ANSI  Std.  239-18 

900.1 n9 


ENERGY  ABSORPTION  FOR  SHORT  DURATION  IMPACTS 


M.  Ziejewski,  B.  Anderson,  M.B.  Rao  and  M.  Hussain 

Technical  Report  96-14 


December  1996 


Center  for  Multivariate  Analysis 
417  Thomas  Building 
Penn  State  University 
University  Park,  PA  16802 


Research  sponsored  by  the  Army  Research  Office  under  Grant  DAAH04-96-1-0082.  The 
United  States  Government  is  authorized  to  reproduce  and  distribute  reprints  for  govern¬ 
mental  purposes  notwithstanding  any  copyright  notation  hereon. 


Typeset  by 


ENERGY  ABSORPTION  FOR  SHORT  DURATION  IMPACTS 

by 

M.  Ziejewski,  N.  Anderson,  M.B.  Rao^  and  M.  Hussain 
North  Dakota  State  University 


ABSTRACT 

The  objective  of  this  study  is  to  relate  energy  absorption  characteristics  to  selected  ma¬ 
terial  properties  and  to  establish  a  methodology  that  allows  one  to  determine  some  of  the 
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Key  Words  and  Phrases:  Energy  absorption.  Density,  Elastic  modulus.  Tangent  modulus. 
Yield  strength,  Poisson  ratio.  Response  surface.  Canonical  analysis. 


^The  research  work  of  M.B.  Rao  is  supported  by  Army  Research  Office  Grant  DAAH04-96-1-0082 


2 


DTIC  QUALTIY  mSPEOIED  9 


Energy  Absorption  for  Short  Duration  Impacts 

Mariusz  Ziejewski,  Beth  Anderson,  Marepaili  Rao,  and  Manzoor  Hussain 

North  Dakota  State  Univ. 


ABSTRACT 

The  objective  of  this  study  is  to  relate  energy  absoiption 
characteristics  to  selected  material  properties  and  to  establish 
a  methodology  that  allows  one  to  determine  some  of  the 
material  properties  for  maximum  energy  absorption.  The 
finite  element  program  DYNA-3D  and  its  associated  pre  and 
post  processors  were  used.  The  model  used  is  a  hollow  square 
column.  Five  properties  of  the  materials  were  included  in 
the  analysis:  (i)  Density  (ii)  Elastic  Modulus  (iii)  Tangent 
Modulus  (hr)  Yield  Strength,  arul  (v)  Poisson  Ratio.  The 
Response  Surface  Method  in  conjunction  with  the  canonical 
analysis  were  employed  to  locate  the  optimum  or  near 
optimum  levels  of  the  properties  and  then  to  determine  the 
equation  of  the  response  surface  in  an  area  near  the  vector  of 
optimum  levels.  For  the  given  levels  of  three  out  of  five 
material  properties  used  in  the  study,  one  can  calculate  the 
remaining  two  material  property  levels  to  achieve  the  near- 
optimal  energy  absoiption. 

INTRODUCTION 

To  design  safer  automobiles,  knowledge  of  the  energy 
absorption  characteristics  of  the  materials  used  in  the 
structural  elements  of  the  automobile  is  very  important. 
There  is  a  substantial  number  of  publications  which  discuss 
the  theory  behind  and  the  experimental  results  of  the 
structural  collapse  of  hollow-thin  wailed  square  tubes 
(McNay,  1988;  Mahmood  et.al.,  1981;  Mahmood  et.al.,  1988, 
Schmueser  et.al.,  1988;  Tundermann,  et.al..  1975).  However, 
there  is  limited  information  on  the  relationship  between  the 
energy  absorption  characteristics  of  materials  and  their 
properties. 

The  general  objective  of  this  study  is  to  relate  energy 
absorption  characteristics  to  selected  material  properties. 

Specifically: 

1 .  To  establish  the  energy  absorption  for  selected  materials. 

2.  To  develop  a  prediction  equation  that  is  capable  of 
determining  energy  absorption  based  on  density,  modulus 


of  elasticity,  tangent  modulus  of  elasticity,  Poisson’s  ratio 
and  yield  strength. 

3.  To  establish  the  level  of  the  selected  material  properties 
for  maximum  energy  absoiption. 

4.  To  establish  a  methodology  that  allows  one  to  determine 
some  of  the  material  properties  for  maximum  energy 
absoiption. 

NUMERICAL  ANALYSIS  USING  DYNA.3D 
For  this  analysis,  the  finite  element  program  DYNA-3D  and 
its  associated  pre  and  post  processors  were  used  on  an  IBM 
RS  6000  Model  970  running  ADC  version  325. 

The  pre  processor  used  for  mesh  generation  was  INGRID 
version  1992a  (Christen  and  Dovey,  1992).  DYNA-3D 
version  3.3.5,  which  is  an  explicit  analysis  code  intended  for 
dynamical  impact  and  crash  analysis  was  utilized  (Whirley 
and  Engelmann,  1993).  The  post  processor  for  the  graphic^ 
output  and  analysis  was  TAURUS  revision  3.0.92  (Spelce. 
1991). 

MODELING  DETAILS 

The  model  represents  a  coluimi,  made  of  1..5  mm  thick  sheet 
metal,  304.8  mm  long  and  having  square  cross-section 
outside  width  dimension  of  69.9  nun.  Load  was  applied  to 
the  upper  50.8  nun  of  the  column  as  a  constant  velocity  field 
which  acts  in  a  direction  parallel  to  the  column's  longitudinal 
axis. 

The  most  unstable  stage  of  the  buckling  is  the  initiation  of 
lateral  deflection.  This  is  numerically  stabilized  in  the  model 
by  using  a  small  crease  or  initial  displacement  in  the  column 
at  the  interface  between  the  upper  and  lower  portions.  This 
crease  starts  the  buckling  in  a  predetermined  direction,  thus 
eliminating  the  initial  numeric  instability. 

The  column  contained  two  planes  of  stmrmetry  along  the  axis 
of  crush.  Therefore  the  crush  behavior  could  be  predicted 
with  a  quarter  model.  A  uniform  mesh  was  used  in  the 
analysis  which  consisted  of  1800  four  node  shell  elements 
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using  three  integration  points  through  the  thickness.  The 
mesh  used  900  elements  for  each  side  of  the  quarter  sector: 
10  elements  for  the  flange  width  and  90  elements  for  the 
flange  length.  Initial  studies  indicated  that  this  high  density 
mesh  was  necessary  to  capture  the  behavior  of  the  comers  and 
hence  the  behavior  of  the  column  as  a  whole. 


The  boundary  conditions  were  specified  in  the  following 
manner.  Symmetry'  boundary  conditions  were  enforced  such 
that  the  quarter  model  might  represent  the  symmetric 
behavior  of  the  full  physical  column.  Note  that  the  column 
beam  walls  fold  onto  themselves  in  a  distorted  sinusoid 
pattern.  To  prevent  the  contacting  surfaces  from  penetrating 
each  other  a  slide  surface  was  defined.  The  particular  slide 
surface  used  was  the  single  surface  contact  (type  4)  slide 
surface.  The  key  feature  of  this  type  of  slide  surface  is  that 
every  node  in  the  definition  is  a  slave  to  all  other  nodes.  The 
advantage  to  using  this  type  of  slide  surface  lies  in  the  fact 
that  any  portion  of  the  defined  area  can  contact  any  other 
portion  without  undesirable  penetration.  The  fixed  end  of  the 
tube  was  constrained  in  all  six  degrees  of  freedom  to  prevent 
rigid  body  motion  when  loaded.  The  moving  end  of  the 
colunm  was  constrained  in  five  degrees  of  freedom  the 
remaining  degree  of  freedom  along  the  longitudinal  axis  of 
the  column  was  left  unconstrained  so  as  to  simulate  the 
moving  head  of  a  hydraulic  press.  All  nodes  on  the  moving 
end,  were  given  a  constant  velocity  of  6934.2  mm/sec  to  crush 
or  collapse  the  column.  The  material  model  used  was  bilinear 
elastic/plastic  with  isotopic  hardening. 


TEST  MATERIALS 

Materials  used  in  the  numerical  analysis  arc  given  in  Table  1. 
Five  properties  of  the  materials  were  used  :  (i)  Density  (ii) 
Elastic  Modulus  (iii)  Tangent  Modulus  (iv)  Yield  Strength, 
and  (v)  Poisson  Ratio. 


The  range  of  the  test  material 
Density  (kg/cm^) 

Elastic  Modulus  (GPa) 
Tangent  Modulus  (MPa) 
Yield  Strength  (MPa) 
Poisson’s  Ratio 


properties  were: 

0.0018  to  0.0091 
44.8  to  244.8 
96.7  to  4468.8 
27.6  to  910.1 
0.264  to  0.35 


STATISTICAL  ANALYSIS  OF  ENERGY  ABSORPTION 
The  main  objective  of  the  statistical  analysis  was  to  locate  the 
optimum  or  near  optimum  levels  of  the  properties  and  then  to 
determine  the  equation  of  the  response  surface  in  an  area  near 
the  vector  of  optimum  levels.  The  Response  Surface  Method 
in  conjunction  with  the  canonical  analysis  was  used  to 
maximize  the  energy  absorption  on  the  basis  of  the  five 
selected  material  properties.  The  Response  Surface  Method 
which  was  first  developed  and  described  by  Box  and  Wilson 
(1951)  is  used  extensively  in  the  field  of  industrial  research. 
Box  and  Wilson  discussed  the  Steepest  Ascent  Method  in  the 
search  for  the  near-stationary  region  surrounding  the 
optimum  and  this  method  was  useful  here  in  searching  for 
optimal  conditions. 


The  Energy  Absorption  was  defined  as  the  response  variable 
and  the  other  five  variables  as  independent  variables  or 
controlled  variables. 


All  the  variables  have  been  scaled  up  or  down  in  order  to 
bring  the  values  of  the  variables  on  par.  The  following  is  a 
description  of  the  scaling  that  has  been  carried  out. 


Variable 

Energy  Absorption 
Density 

Elastic  Modulus 
Tangent  Modulus 
Yield  Strength 
Poisson  Ratio 


Scaling  Factor 
1/10,000 
1/1,000 

1/1,000,000.000 

1/10,000,000 

1/1,000,000 
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In  order  to  determine  the  equation  of  the  response  surface,  a 
simple  model  based  on  five  independent  variables  was  defined 
as: 


Y  =  Bo+B,  X,  +B2X2  +  B,X,  +  B4X,  +  B5X5  +  e,  (1) 


where 

Y  =  Energy  Absorption  is  measured  in  J, 

Xi  =  Density  of  the  material  is  measured  in  kg/cm’. 

X2  ==  Elastic  Modulus  is  measured  in  Pa, 

Xb  =  Tangent  Modulus  is  measured  in  Pa, 

X4  =  Yield  Strength  is  measured  in  Pa 
Xs  =  Poisson  Ratio. 

It  was  assumed  that  e  is  the  random  error  having  mean  zero 
and  variance  .  The  estimates  of  the  B’s  are  to  be 
determined  by  the  method  of  least  squares  which  minimizes 
the  sum  of  squares  of  the  errors  e. 

Since  the  preliminary  analysis  indicated  that  the  surface  was 
not  planar,  a  higher-order  equation  in  the  variables  was  a 
more  suitable  model  to  find  the  optimal  Ie\'els  of  properties. 
The  model  was  defined  as: 

Y  =  Bo+  B,  X,  +  B2  Xj  +  B3  Xj  +  B4  X4  +  Bs  Xs  +  B„  X,’  + 
B22  X2^+  B33  Xj^+  B44  X4^  +  B55  Xs^  +  B,2  X,  X2  +  Bi3  X,  X3 
+  Bu  X,  X4  +  B,5  X  iXs  +  B23  X2  X3  +  B24  X2X4  +  B25  X2  Xs 

+  B34X3X4  +  B35X3X5  +  B45X4X5  +  e.  (2) 

The  time  determinations  are  done  using  the  canonical 
analysis  method. 

The  canonical  method  is  a  method  of  rewriting  a  fitted  second 
degree  equation  by  changing  the  origin  such  that  first-order 
terms  and  cross-product  terms  have  been  removed.  It  is  also 
called  B  canonical  form.  Box  and  Draper  (1987)  have  given 
full  details  of  this  procedure.  If  the  determined  stationary 
point  is  a  saddle  point  it  allows  considerable  choice  of 
operating  conditions  (Box,  1954). 
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RESULTS 

Typical  partial  output  from  a  compression  test  is  given  in 
Figures  1  and  2. 


t  =•  0  msec  t  “  30  msec  t  ”  60  msec 

<l  =  0mm  d=125.7mm  d  =  231.9mm 


Figure  1.  Deformed  geometries  of  column  during 
compression. 


Deformation  (mm) 

Figure  2.  Force-deformation  curve  of  typical  coluiim 
compression. 

Using  the  regression  model  given  above  (1)  the  results  are: 

Y  =-21.2  -0.5  X,  +0.1X2  + 

0.07X3+0.07X4  +  5.6X5  (3) 

with  F  =  691.0  and  conesponding  p=0.0001  and  R^=0.99. 

The  Response  Surface  Model  (2)  for  the  given  data  is: 

Y=  8.3  +  61.9  X,  -  2.5  X2  +  1.4  Xj  -  0.2  X4  -  9.8  X5+  3.5 
X,-  -  0.5  X,  X2  +  0.01  X2-+O.4  XiX,  -0.01  X2X3-O.OOO4  X,^  - 
0.03  X,X4  +  0.001  X2X4  -  0.00002  X3X4  -  0.00006  X4--  16.9 
X,Xs  +  0.7  X2X5  -  0.4  X3X5  +  0.09  X4X5  +  1.8  X5- ,  (4) 

withR-=  1.0. 


The  stationary  point  T,  x,  =  (x,,  .x,,  X3,  X4.  X5).  The  fitted 
model  (4)  can  be  written  in  the  following  way: 


x,' 

X2 

x  = 

X3 

(5) 

X4 

X5. 

B,' 

6L9‘ 

B2 

-2.5 

B  = 

B3 

= 

1.4 

f 

B4 

-0.2 

Bs. 

-9.8 

Bn 

-Bi2 

2 

-  Bn 
2 

-Bi4 

2 

-Bi5 

2 

}b,2 

B22  “  B23 

-B24 

2 

-B25 

2 

H  = 

jB,3 

J  B23  B33 

-B34 

2 

-B35 

2 

-Bi4 

2 

jB24  jBs4 

B44 

^B43 

2 

-B,j 

.2 

-B25 

2 

-B45 

2 

B55 

■  3.5 

-0.2 

0.2 

-0.01 

-8.5 

-0.2 

0.01 

-0.007 

0.0007 

0.4 

0.2 

-0.007 

-0.0004 

-0.00001 

-0.2 

-0.01 

0.0007  -0.00001 

-0.00006 

0.05 

-8.5 

0.4 

-0.2 

0.05 

1.8 

Note  that  if  all  the  eigenvalues  of  H  are  positive,  Y  will  be 
minimum  at  the  stationary  point.  If  all  the  eigenvalues  are 

negative,  Y  will  be  maximum  at  the  stationary  point.  If  the 
eigenvalues  are  of  mixed  signs,  the  stationary  point  is  a 
saddle  point. 

Note  that  ,  i=l,2,...,5  are  the  eigenvalues  and  m„ 
i=I,2,...,5  are  the  corresponding  eigenvectors  of  the  matrix  H. 
For  the  given  problem,  the  eigenvalues  in  order  of  absolute 
magnitude,  and  their  associated  eigenveaors  are  as  follows: 

X  1=11.2,  mi=(0.7,  -0.04,0.02,-0.004,  -0.7), 

X2  =  -5.8,  m2=(  0.71, -0.02,  0.003, -0.004.  0.7), 

X  3  = -0.007,  m3=(-0.03, -0.3,  0.9, -0.1,  0.02), 

X  4  =  0.0003,  m4=(0.03,  0.9,  0.3, -0.3,  -0.009), 

X5  =  -0.0,  m5=(0.01,  0.3,  0.2.  0.9,  0.00008). 

We  \vtite 
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X 


Xi 

X2 

X3 

X4 

.X5 


=  M(X-X,), 


where 


There  is  considerable  flexibility  in  the  levels  of  material 
properties  which  can  give  rise  to  the  value  Y The  model 

Y  can  be  written  in  a  new  coordinate  system  as  follows: 

Y  =  Yj  +X,X?+....+>.5X^ 


Y=14366+1125?  -5S5^ 

-00075^ +0.00033^-05^ 


m, 

1112 


M  = 


mj 


m4 


mj 


X  is  defined  in  (S)  and  X,  will  be  evaluated  as 
-2HXs  =  B 


i.e.. 


X,  =--h'’b  = 


263.4 

6212.5 
5037.7 

20392.5 
5.8 


The  dominant  eigenvalues  are  X]  and  X.2.  The  remaining 
eigenvalues  can  be  ignored.  The  model  can  be  rewritten  as 

YsY,+11.2X,^- 5.83(2^  (12) 

Any  levels  of  Xi,  X2,  X3,  X4,  and  X5  satisfying 

3C,  =0  and  3(2  =0. 

would  give  the  near*optimal  value  Y Depending  on  the 
given  circumstances,  the  engineers  have  a  good  degree  of 
flexibility  in  choosing  the  levels  of  X’s. 

ENGINEERING  APPLICATION  OF  RESULTS 

There  are  two  ways  how  the  results  from  this  analysis  can  be 

used  by  a  design  engineer; 


which  is  called  the  coordinates  of  the  stationary  point. 
This  gives, 

3C,  =  ().7X,  -0.04X2  +0.02Xj 
-0.004X4-0.7X5+2.4  (6) 

3(2  =  0.6X,  -0.02X2  +0.003Xj 
-0.004X4+0.7X5-2.94  (7) 

Xj  =  -0.03X,  -0.3Xj  +0.9X3 
-0.1X4+0.02X5-6.4  (8) 

3(4  =  0.03X,  +0.9X2  +0.3X3 
-0.3X4  -0.009X5  +20.7  (9) 

3(5  =  o.oix,  +0.3X2  +0.2X3 
+0.9X4+0.00008X5-21906.6  (10) 

The  value  of  Y  at  the  stationary  point  is  given  by 

Y,  =B„ +4  X,B  =  1436.6. 


1.  Using  equation  (3)  or  (4)  one  can  determine  the 
estimated  energy  absorption  with  the  knowledge  of  the 
selected  material  properties  in  this  study. 

2.  Setting  the  equations  (6)  and  (7)  to  zero  and  setting  three 
out  of  five  material  properties  at  desired  levels,  one  can 
calculate  the  remaining  two  variables  for  near-optimal 
energy  absorption. 

For  example,  for  fixed  values  of  density,  elastic  modulus,  and 
tangent  modulus,  the  information  on  what  value  of  the  yield 
strength  and  Poisson's  ratio  required  to  obtain  the  near- 
optimal  energy  absorption  can  be  obtained  from  calculations 
using  equations  (6)  and  (7).  If  the  obtained  values  for  the 
yield  strength  or  Poisson’s  ratio  are  not  achievable  from  a 
practical  point  of  view,  one  can  change  one  of  the  previously 
fixed  properties  and  see  what  effect  it  would  have  on  the 
values  of  the  calculated  properties. 

The  analysis  of  equation  (12)  shows  that  there  is  also  a 
possibility  of  higher  energy  absorption  than  calculated  by 
setting  equations  (6)  and  (7)  to  zero.  However,  the  obtained 
results  are  beyond  the  experimental  range  used  in  this  study. 

SUMMARY  AND  CONCLUSIONS 
1.  Within  the  tested  experimental  range  of  properties,  the 
developed  regression  model  and  response  surface  model 
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can  be  used  to  determine  the  energy  absorption 
characteristics. 

2.  For  the  given  levels  of  three  out  of  five  material 
properties  used  in  the  study,  one  can  calculate  the 
remaining  two  material  property  levels  to  achieve  the 
near-optimal  energy  absorption. 

3.  Using  the  developed  technique  one  can  select  the  easiest 
achievable  combination  of  the  chosen  material  properties. 

4.  There  is  considerable  flexibility  in  the  levels  of  the  tested 
material  properties  which  can  assure  near-optimal  energy 
absorption. 
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